^READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMOrR 


MONITORING  AGENCY  NAME  &  AOORESSfll  diffefoni  from  Controlling  Olfice}  IS.  SECURITY  CLASS,  fof  thie  feportj 

UNCLASSIFIED 


ISa.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


17.  DISTRIBUTION  STATEMENT  (of  tho  abafracf  onterod  in  Block  20,  If  different  from  Report) 


t9.  KEY  WORDS  (Continue  on  reveree  side  If  necessery  end  identify  by  block  number) 


2^.  abstract  (Continue  on  reverse  side  If  necesmmry  and  Identify  by  block  number) 

""^Output  feedback  with  a  series  compensator  is  considered  for 
controlling  SISO  discrete  systems.  MSE  criterion  of  design  is 
used  in  tracking  a  given  input  sequence.  The  plant  can  be 
unstable  and/or  nonminimum  phase,  while  the  compensator  is 
constrained  to  be  physically  realizable  and  closed-loop  configura¬ 
tion  is  both  stable  and  physically  realizable.  .Using  variational 
methods  the  required  compensator  is  found,  and  shown  to  be  the 
i^ecessary  and  sufficient  optimal  solution.  4 — _ 

bD7S-,.  1473  ^(9  3.  _  UNCLASSIFIED  _ ' 

security  cl  ASSiFfCATjON  OF  ThIS  PAGE  /Wien  DaJa  Entered) 


AFOSR-TR-  7  9  -  1284 


DESIGN  OF  OPTIMAL  TRACKING  DISCRETE 
CONTROL  SYSTEM^ 


BY 


Y.  Bar-Ness*  j 

Lefschetz  Center  for  Dynamical  Systems'^ 
Division  of  Applied  Mathematics 
Brown  University 
PROVIDENCE/  Rhode  Island  02912 


AND 


J.  Feinstein 
School  of  Engineering 
Tel  Aviv  University 
ISRAEL 


JuLY/  1979 


This  research  was  supported  by  the  Air  Force  Office  of 
Scientific  Research  under  #-AF0SR  77-3063. 

Author  on  sabbatical  leave  from  the  School  of  Engineering, 

Tel  Aviv  University,  Tel  Aviv,  Israel. 

79  12  18  018 

J»roved  for  public  release; 

distribution  unlimited. 


DESIGN  OF  OPTIMAL  TRACKING  DISCRETE 


CONTROL  SYSTEM 

Y.  Bar-Ness  and  J.  Feinstein 


Abstract 

Output  feedback  with  a  series  compensator  is  considered 
for  controlling  SISO  discrete  systems,  MSE  criterion  of 
design  is  used  in  tracking  a  given  input  sequence.  The  plant 
can  be  unstable  and/or  nonminimum  phase,  while  the  compensator 
is  constrained  to  be  physically  realizable  and  closed-loop 
configuration  is  both  stable  and  physically  realizable.  Using 
variational  methods  the  required  compensator  is  found,  and 
shown  to  be  the  necessary  and  sufficient  optimal  solution. 
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I .  Introduction . 

The  design  of  a  discrete  system  to  track  a  given  input 
sequence  is  an  important  problem.  The  use  of  the  minimum 
squared  error  (MSE)  as  a  criterion  of  performance  in  such 
problems  is  well  accepted.  Kue^^^  in  treating  this  problem, 
did  not  take  into  consideration  either  the  system  configuration 
or  the  plant  properties,  so  that  in  his  solution  the  closed- 
loop  may  be  unstable  and  the  compensator  obtained  might  not  be 
physically  realizable. 

In  this  paper,  output  feedback  together  with  a  series 
compensator  is  considered  for  controlling  the  system  to  obtain 
the  desired  objective.  This  procedure  is  particularly  suitable 
if  the  required  compensator  is  constrained  to  be  physically 
realizable^  ^  and  the  resultant  closed-loop  system  is  to  be 
both  stable  and  physically  realizable. 

Consider  the  physically  realizable  plant  transfer  function 


G(z)  = 

q'^(z)q’  (z) 


where  p^  and  q^  are  polynomials  whose  zeros  are  located 
inside  the  unit  circle  and  p  and  q  are  polynomials  whose 
zeros  are  located  outside  this  circle. 


System  is  physically  realizable  iff  the  difference  between 
the  degrees  of  its  transfer  function's  denominator  and 
nominator,  ^  >  is  not  less  than  zero. 
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For  the  given  input  sequence  U  and  its  corresponding 

k 

output  sequence  Wj^  the  squared  error  is  given  by, 


J  =  I  el  =  I  (U.-W  (2) 

k=0  ^  k=0  ^  ^ 


Using  Parseval's  theorem,  we  get 


J 


1 

2nj 


* 

0  E(z)E(z 


- 1-.  dz 
^  z 


(3) 


where 


ECz)  =  UCz)  -  W(z)  =  I  e  z'^  (4) 

k=0  ^ 

is  the  z  transform  of  the  error  sequence. 

Denoting  the  closed-loop  transfer  function  by  T(z),  we 

have 


E(z)  =  (l-T(z))U(z) 


C5) 


with 


TCz) 


G(z)D(z) 

l+G(z)DCz) 


(6) 


where  DCz)  is  the  compensator  transfer  function  that  has  to 
be  obtained  such  that  the  squared  error  is  minimized. 
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To  solve  the  minimization  problem,  it  is  possible  to  use 
the  Wiener  filter  approach  and  find  TCz).  T(z)  may  even  be 
required  to  be  a  physically  realizable  solution.  Nevertheless, 
with  T(z)  available,  the  required  compensator  is  obtained 
from  equation  (6) 


D(z) 


=  1  TCz) 

Gl^  1-TCz)  • 


(7) 


Clearly,  DCz)  must  cancel  some  of  the  zeros  or  poles  of  the 
plant  GCz),  which  might  be  outside  the  unit  circle,  thus 
yielding  an  unstable  closed-loop  system.  To  cope  with  the 
problem  of  zero  cancellation  in  the  continuous  case, 
Newton^^^  and,  later,  Bongiorno^^^  and  Weston^^^^  restricted 
the  plant  to  be  asymptotically  stable.  Instead,  as  is  done 
in  this  paper,  it  is  reasonable  to  require  that  D(z)  is 
allowed  to  cancel  only  those  poles  and  zeros  of  G(z)  which 
are  inside  the  unit  circle.  That  is 

D(z)  =  .  (8) 

P  (z)B(z) 


where  “(z)  and  B(z)  are  polynomials,  yet  to  be  determined. 

To  obtain  these  polynomials  we  define  an  appropriate 
variational  problem  on  a  space  of  polynomials  such  that  it 
takes  into  account  the  requirements , that  DCz)  be  physically 
realizable  and  that  TCz)  is  both  stable  and  physically 
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realizable.  From  the  solution  of  the  variation  problem  we 
obtain  the  necessary  and  sufficient  condition  which  D(z)  - 
in  particular  and  3(z)  -  has  to  satisfy.  It  turns 

out  that  for  D(z)  to  be  the  required  optimal  solution,  »(z) 
and  3Cz)  must  satisfy  some  polynomial  equation  whose  para¬ 
meters  are  easily  obtained  from  the  input  transform  UCz) 
and  the  plant  G(z)  -  see  Theorem  4.6.  Furthermore,  we  show 
that  under  a  certain  polynomial  degree  condition,  the  solution 
of  this  polynomial  equation,  “(z)  and  B(z) ,  is  unique  -  see 
Lemma  4.4.  With  this,  D(z)  is  the  required  optimal  compensator 
see  Theorem  4.7.  Finally,  we  calculate  the  value  of  the 
optimal  performance  for  comparison  purpose.  Two  examples 
supplement  the  theoretical  results. 

Although  we  consider  the  case  of  a  given  input  sequence, 
nevertheless,  it  is  possible  to  extend  this  work  to  the  case 
when  the  input  includes  noise  or  distrubances  with  a  given 
spectral  density. 

Youla  et  al^^^  treated  a  related  problem  for  the  con¬ 
tinuous  case,  and  (as  in  this  work)  the  plant  was  allowed  to 
be  unstable  and/or  nominimum  phase.  In  addition  to  the  funda¬ 
mental  difference  between  the  discrete  and  continuous  systems, 
Youla' s  work  differs  from  ours  in  several  respects.  The 
approaches  are  quite  different.  Youla's  paper  uses  spectral 
factorization  and  some  of  the  implicit  assumption,  might  be 
rather  restrictive  (e . g .  assumption  6).  Secondly  the  resultant 
controller  might  be  physically  unrealizable  depending  on  the 


given  input  or  plant,  a  situation  which  might  also  be  con¬ 
sidered  restrictive,  particularly  if  flexible  implimentation 
of  this  controller  is  required  (see  [5]  foortnote  on  page  7). 
Finally,  Theorem  1  of  Youla  constitutes  only  a  necessary 
condition.  Another  reference  which  may  be  mentioned  is 
Volgin^^^,  who  discussed  the  problem  without  adequately 
solving  it. 

II.  Stating  the  Problem. 

Consider  the  discrete  control  system  with  output  feedback 
and  a  series  compensator.  Let  the  plant  G(z)  be  given  by 
equation  (1).  For  a  given  input  sequence  Uj^,  it  is  required 
to  find  the  polynomials  “(z)  and  3(z)  of  the  compensator 
DCz)  (equation  (8)),  such  that: 

(1)  D(z)  is  physically  realizable. 

(2)  The  closed-loop  transfer  function  T(z)  (equation  (6)) 
is  stable  and  physicall  realizable. 

(3)  The  sequared  error  between  the  given  input  sequence 

Uj^  and  its  corresponding  output  sequence  is 

minimized.  (Equivalently,  the  cost  functional  of 
equation  (3]  is  minimized.) 

III.  On  the  Constraints  of  the  Variational  Problem. 

First,  we  consider  some  results  related  to  the  constraints 
of  the  variational  problem.  In  terms  of  the  polynomials  cx(z) 
and  3(z) ,  the  closed-loop  transfer  function  is  given  by 


T(z)  = 


Cz)cx(z) 


where , 

v(z)  =  p'(z)ot(z)  +q  Cz)3(z),  (10) 

Propcsition  3.1:  For  T(z)  to  be  physically  realizable  it  is 
necessary  and  sufficient  that  the  degree  6^  of  the  polynomial 
v(z)  in  equation  (10)  satisfies 


6  =  max  (6  _  ,  <5  _  ) 

p'a  q'3 


Proof:  From  equation  (9),  if  T(z)  is  physically  realizable, 


then 


£  A  6  -  6  >0 

^  ^  pa 


^  (S  >  6 

V  - 

pa 


From  equation  (10) 


6  =  max  (6  _,6  )  if  6_ 

paqu  pa  qB 


<6_=6_  if6_=6_ 

p'a  q  B  p'a  q"B 


and  together  with  equation  (12)  the  result  follows.  The  other 


r 


direction  of  our  claim  is  obvious. 

Proposition  3.2:  Assume  that  the  plant  GCz)  is  physically 
realizable,  then 

(a)  For  both  the  closed-loop  transfer  function  T(z)  and 
the  compensator  Dfs)  to  be  physically  realizable  it  is 
necessary  that 


in  which  case, 


6=6  >6 
^  q'6  q'“ 


(13) 


I 


T 


(14) 


where 


is  the  excess  denominator- nomina tor  zeros. 


(b)  On  the  other  hand,  if 


6  =  6  _  >  6  _ 
q  3  ’  p 


(13) 


and 


>  Hq. 


(15) 


Then  both  T(z)  and  D(z)  are  physicall  realizable. 
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Proof :  (a)  Using  the  fact  that  GCs)  is  physically  realizable 

we  have  from  equation  (1) 

Jl„A6^+6  -(6^  +  6  )>0  (16) 

q  q'  P  P 

Similarly  D(z)  is  physically  realizable  implies 

*n  A  «  *  *  ■  (4  0-  (17) 

p  q 

These  two  equations  together  yield 

6  .  +  6g  =  6  .  >  6  _  =  6  _  +  6^.  (18) 

q  ^  q  8  p  Cl  p 

But,  by  Proposition  3.1,  for  T(z)  to  be  physically  realizable 
it  is  necessary  that 


6  =  max(<5  6  )  Cf9) 

p'oi  q’8 

and,  together  with  equation  CIS),  the  claim  follows. 

Furthermore,  adding  equations  (16)  and  (17),  we  have 


6 

q*  8 


6 

p'a 


6  - 
V 


6  A  I 
p’a  ^ 


(b)  From  equation  (13)  and  Proposition  3.1,  T(z)  is 


physically  realizable.  Also,  by  equation  (13) 


and  together  with  equation  (15),  using  equation  (16)  for  , 

U 

yield 


«  .  -  <50  -  (6  .  >  6  ^  .  6  _  -  (6  ,  .  6  J  >  0 

q  p  q  q  p  p 


^  S  ‘  ^  >  0. 

p  q 


Hence,  using  equation  (8),  D(z)  is  physically  realizable. 


IV .  The  Solution  of  the  Variational  Problem 

For  any  given  polynomial  x(z)  of  degree  6^,  let 

6 

x(z)  A  z  ^  x(z  ■^)  .  (20) 


n 

If  x(z)  =  I  a.z^,  a  ^0  6  =  n,  then  clearly 

i=0  ^  n  X 


11  ,  11 

x(z)  =  z  I  a.z^  =  I  a.z^  Also,  if  x(z)  has  no 

i=0  ^  i=0  ^ 


at  z  =  0  (i.e.,  a^  ^  0),  then 


zeros 


6~  = 
X 


6  . 

X 


(21) 


Proposition  4.1:  Let  $(z)  be  a  rational  function  of  z  whose 
aero  and  poles  are  symmetrically  located  with  respect  to  the 


unity  circle  ,  i . e .  , 


<{>(z)  =  UCz)UCz‘-^)  . 


Then,  it  is  possible  to  factor  iI>Cz) 


(z)'^  Cz) 
p'^Cz)P  (z) 


where  the  zeros  of  and  P^  are  located  inside  the  unit 

circle  and  those  of  t  and  P  are  outside  it,  such  that 


X  (z)  =  t'(z) 
^■"(z)  =  P'Cz) 


6=6 
+  + 

T  P 


Proof :  See  Appendix  A,  claim  (A-1). 

We  are  nov;  in  a  position  to  set  the  proper  variational 
problem.  Substituting  equation  (5)  into  equation  (3)  yields 


J  =  ^  ^  (l-T(z))(l-T(z'^))'I>(z)^ 


where  T(z)  is  given  by  equations  C9)  and  (10)  and  $(z) 


by  equation  (22).  Therefore, 
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6  > 
a  - 


(30) 


and 


6  =  6  _  >  . 

^  q  3  p  “ 


(31) 


then  for  (“ , 3)  £  P  to  be  extremum  of  the  cost  functional  of 
equation  (26) ,  it  is  necessary  that 


c(z)  =  q'(z)3(z)p'Cz)q'C2)'i:'^(z)T'(2)z 
v(z) [v(z)f p^(z)p'  (z)z 


(32) 


where 


v(2)  =  p  (z)a(z)  +  q  (z)3(z)  (33) 

has  no  poles  inside  the  unit  circle. 

Proof :  For  J  to  have  an  extremum  at  (*^ ,  3)  ^  P,  the 

variation  of  equation  (29)  must  vanish  for  every  polynomial 
increment  hj(z)  and  h2(z).  Define 

I 

A(z)  A  p  (z)q  (z)  3(z)t‘*'(z)t' (z)p  (z)z  ^/d(z) 


and 
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B(z)  A  zv(z)(}>^(z)/dCz) 

so  that  ACz)  and  B(z)  have  no  common  zeros  (i.e.,  d(z)  is 
the  corresponding  g.c.d.)-  Also,  let 

p(z)  A  C3(z)h^Cz)  '  ct(z)h2(z)). 

The  polynomials  v(z)  and  p’Cz)  have  no  zeros  inside  the 
unit  circle  and  can  be  disregarded.  Clearly, 

5  1  *  ‘v  ” 

For  a  given  ^Cz)  and  3Cz)  ,  h^Cz)  and  ^2^^^  selected 

such  that  M(z)  could  have  zero  of  order  ^  5  anywhere 

in  the  complex  plane,  and  particularly,  inside  the  unit  circle. 

From  equations  (30),  (31),  (34)  and  since  <5  ^.  =  <5  ^ ,  by 

P 

equation  (25)  we  deduce  that  Lemma  B-2 

our  claim  follows. 

Lemma  4.3:  Assume  that  the  closed-loop  transfer  function  T(z) 
is  stable;  then  for  every  >  0,  for  c(z)  to  be  analytic 
inside  the  unit  circle,  it  is  necessary  that 

I 

v(z)  =  t'^(z)p' (z)q' (z)z  ^  (35) 

6(z)  =  zp'^(z)Tr(z) 


(36) 


where  Tt(z)  is  some  polynomial  which  is  together  with  a(z) 
must  satisfy  the  following  equation 


p(z)a(z)  +  q‘ Cz)  •zP’^(z)Tr(z)  =  t'^(z)p' Cz)q' (z)  z  (37) 

If,  furthermore. 


6=6  >6 

q'$  ■  p‘“ 


then  it  is  necessary  that 


(38) 


and 


(39) 


6^  =  6  ^  +  6  _  (40) 

T  q 

where  A  6  -  6 

q'6  p'a 

On  the  other  hand,  if  “(z)  and  ''^(z)  are  a  polynomial 
pair  solution  to  equation  (37),  then  with  3(z)  as  in  (36), 
c(z)  is  analytic  inside  the  unit  circle. 

Proof :  For  T(z)  to  be  stable  v(z)  must  contain  only  zeros 

inside  the  unit  circle.  Since  Ij  might  be  zero  then  for 
c(z)  to  be  analytic  inside  the  unit  circle  we  must  have 
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v(z)  zp‘^(z)'fr  (z)  =  3(z)p  (z)q  (z)i‘^(z)z  ^  (41) 


for  some  polynomial  ffCz).  Substituting  equation  (33),  and  re¬ 
arranging  terms,  we  have 


3(z)  [p' (z)q' (z)t'^(z)z  ^  -  zq  (z)  (z)  tt  (z)  ] 


=  p  (z)u(z)if(z)P  (z)z 


(42) 


Since  3(z)  and  p  (z)ci(z)  are  relatively  prime  and  ''^(z)  is 
arbitrary,  equation  (36)  follows.  And  consequently,  from 
equation  (42)  ''t(z)  must  satisfy  equation  (37).  Substituting 

equation  (36)  in  equation  (37)  results  in  equation  (35). 

Furthermore,  using  equation  (38)  we  have 


(43) 


Since  p*(z)  and  q  (z)  have  no  zeros  inside  the  unit  circle, 
6^  =  ^  .  and  6^  ~  ^  and  from  equation  (35)  we  have 

p’  p"  q*  q' 


T  p  q 


(44) 


Also,  from  equation  (36)  we  have 


6=1.6.  6,. 

(f 


(45) 
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Combining  equations  (43),  (44)  and  (45)  and  using  equation  (25), 
equation  (39)  follows.  Also,  from  equation  (9) 


A  6 

T  —  V 


6 

p'a 


f 


and  substituting  for  from  equation  (44) ,  equation  (40) 

follows  after  simple  manipulations. 

On  the  otherhand,  if  we  pick  ^^(z)  and  “(z)  to  be  a 
solution  of  equation  (37)  and  take  3(z)  as  in  equation  (36) , 
then  clearly  v(z)  will  have  the  form  of  equation  (35)  and 
using  equation  (32)  c(z)  is  obviously  analytic  inside  the 
unit  circle. 


Lemma  4.4;  The  solution  pair  of  polynomials  '''(z)  and  ci(z) 
of  equation  (37) , 

£ 

p'(z)a(z)  +  q' (z)zp'^(z)TT(z)  =  t'^(z)p*  (z)q‘ (z)  ^  (37) 

which  satisfy 


6,  =  6  .  -H  6  -  1 

P 


^cx  =  ^  ^  - 

q 


is  unique  for  all  >  6^ 


(39) 
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+  T  + 

Proof:  Assume  that  t  Cz)p  Cz)q  (z)z  and  zp  (z)  are 

relatively  prime.  To  solve  equation  (37)  for  the  unknown 
polynomials  o^Cz)  and  ’’^Cz),  notice  that  the  number  of 
unknowns  (the  coefficients  of  these  polynomials)  is 
ni  =  <5^  +  +  2 .  Using  equations  (39)  and  (40) 

n^=6_+6^+6^+6_+l.  (41) 

P  q 

On  the  otherhand,  the  number  of  scalar  equations  obtained  by 
equating  the  coefficients  of  different  powers  is  given  by  the 
degree  of  the  polynomial  on  the  left 

n-  =  6  ^  +  6  +6  +6+1.  (42) 

^  P  ?' 

Since,  6^  =  6  and  6=6  it  follows  that  n,  =  n-. 

p'  p'  q'  q' 

Now,  the  polynomials  p  (z)  and  q  (z)zP  (z)  are 

relatively  prime  and  hence  their  resultant  is  nonzero,  and  the 

matrix  of  the  coefficients  is  nonsingular.  his  ensures  the 

existence  and  uniqueness  of  the  solution. 

+  --  + 

If  the  polynomials  t  (z)p  (z)q  (z)z  and  zP  (z)  are 

not  relatively  prime,  we  denote  their  g.c.d.  by  d(z).  Let 

a(2)  =  d(z)a'(z) 

z  P  (z)  =  d(z)  P  (z) .  (43) 


=  dCz)v'  (z)  . 


(44) 


■[■^(z)!!' (z)q' (z)z  ^ 


Equation  (37)  thus  becomes 

p  (z)a'(z)  +  q  '  (z)p  '  (z)  TT  (z)  =  v' (z)  , 

This  leaves  the  number  of  equations  and  unknowns  equal  and, 
therefore,  the  previous  argument  of  existence  and  uniqueness 
for  the  pair  “(z)  and  ''^(z)  follows. 


Lemma  4.5:  If 


c(z)  =  q'(z)3(z)p'(z)q'(z)t'^(z)t'(z)z 
v(z)  [v(t)  ]^p'*' (z)p‘  (z)z 


(32) 


is  analytic  inside  the  unit  circle,  then  the  second  variation 
of  the  cost  functional  of  equation  (26)  is  given  by 


y2j[h,,h2]  =  ^ 

^  ^  J  J  [v(z)]^  [v(z)]^ 

•  [g(z)h^(z)  -  <5;(z)h2(z)]  [6(z)h^(z) 


(45) 

“(z)h2(z)] 


and  is  positive  semidefinite. 

Proof :  Notice  that  the  integrand  of  the  first  integral  of 

equation  (B-ll)  in  Appendix  B  equals  c(z)/v(z).  Therefore,  it 
is  analytic  inside  the  unit  circle  and  the  integral  is  zero  for 


k 
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every  selection  of  h^^Cz)  and  h2(2).  The  second  integral  of 
equation  (B-11)  which  is  given  in  equation  (45)  can  be  rewritten 
in  the  form  (see  also  equation  (B-5)) 


u2j[h^.h2]  =  X  ^  u[T(z)]p[T(z'^)]<I>(z)- 


Using  Parseval's  theorem,  equation  (46)  can  be  represented  as 

2 

a  sum  of  squares,  hence  M  J[  ]  is  positive  semidef inite . 
Finally,  our  results  can  be  summarized  as  follows. 


Theorem  4.6  (necessary  condition):  Let  the  input  sequence  for 
the  discrete  control  system  with  output  feedback  be  Uj^  and  its 
z-transform  U(z).  Define 


$(z)  A  U(z)’j'^(z)  = 

P  (z)P  (z) 


If  the  plant  is  physically  realizable  and  given  by 


(z)p' (z) 


q  (z)q'(z) 


then  for  the  compensator 


D(z)  =  ^ 

p  (z)B(z) 


to  be  the  optimal  solution  to  the  problem  stated  in  Section  II, 
it  is  necessary  that 
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BCz)  =  zP  (z)^Cz) 


+  -^X 

v(z)  =  T  (z)p  (z)q  Cz)z 


(47) 

(48) 


where 


v(z)  =  p'(z)a(z)  +  q  (z)6(z), 


(49) 


and  7r(z)  is  a  solution  to  the  polynomial  equation 


p‘(z)<Jt(z)  +  q'(z)zp'^  (z)Ti(z)  =  (z)p' (z)q' (z)  z 


with 


6^  =  6  .  +  il  -  1,  (£  A  6  ,  -  6  .  ) 

p  *  q  3  p  a 


6,  =  6  +6 
a  + 

^  q 


and  where  q^,q  ,p^  and  p  are  defined  by  equation  (20) 


Proof :  Using  Proposition  3.2,  we  have 


6v  =  6  >  <5  _ 

q  3  p  “ 


and  using  Lemma  4.2  and  4.3  our  claim  follows. 


V 


•■F  I"*' 
A 
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Theorem  4 . 7  (sufficient  condition):  For  the  discrete  control 
system  with  output  feedback  and  physically  realizable  plant 
as  in  Theorem  4.6,  let  ifCa)  and  a(z)  be  the  solution  of 

I 

p'(z)ct(z)  +  q  (z)zp'^(z)7r(z)  =  t‘^(z)p  (z)q  (z)z  ^  (37) 

satisfying 


6 


n 


+ 


1 


6  = 
a 


6 

q 


where  >  0  and  the  different  polynomials  are  specified  in 
Theorem  3.9, 

Let 


B(Z)  =  Zp'^(z)TT(z) 


(50) 


then  using 


D(z)  .  (51) 

P  (z)B(z) 

as  a  series  compensator,  the  closed-loop  system  is  stable,  the 
first  variation  will  vanish  and  the  second  variation  is  positive 
semi-definite  for  all  increments  in  “(z)  and  B(z) .  If, 
furthermore,  is  chosen  so  that  8.^  >  then  D(z)  is  also 
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physically  realizable. 


Proof :  First  using  Lemma  4.4,  the  pair  and  ''^(z)  is 

unique.  From  equation  (8) 


T(z) 


^  p  (z)aCz) 
v(t) 


where  v(z)  is 

£ 

v(z)  =  'c’^(z)p  (z)q  (z)z  ^  (52) 

due  to  the  way  d(z)  is  selected.  Hence,  T(z)  is  stable,  and 

6  =  6^  +  6  +6  +Jt„. 

V  +  -  -  T 

q  r 

By  Lemma  4.3,  c(z)  (equation  (32))  is  analytic  inside  the  in 
unit  circle  and,  by  equation  (29),  the  first  variation  is  zero. 
Also,  by  Lemma  4.5  the  second  variation  of  equation  (53)  is 
positive  semidefinitr  for  all  increments  of  the  polynomials  0‘(z) 
and  B(z) , 

Now,  by  our  choice  of  7T(z),a(z)  and  consequently  of 

B(z) 


6  _  =  6  _  +  1  +  6  ^  +  (6  _  +  £  -  1) 
q  q  p  p 


^  =6  +6^+6 
0  -  -  + 
pap  X  q 


Therefore,  using  ^  ^  +  and  £„  >  0  we  have 

6  =  6  _  >  S  .  . 

q  3  p  ot 

and  if  we  choose  2.^  >  then,  by  Proposition  3.2  the 

compensator  D(z)  is  physically  realizable. 

V .  The  Minimum  Square  Error. 

For  the  unconstrained  problem  we  could  have  used  the 
Wiener  solution  to  obtain  the  open-loop  otpimal  filter.  For 
the  case  of  simple  poles,  the  optimal  filter  transfer  function 
(our  closed-loop  transfer  function)  obtained  using  partial 
fraction  expansion^^^  has  no  excess  of  poles  over  zeros,  i.e., 

=  0. 

Realizing  this  optimal  Wiener  filter  T(z)  using  series 
compensator  configuration,  we  have  that 

=  _ 3: _  -T I- 

G(z)  l-T(z)  • 

Hence,  if  >  0 ,  2,^  =  and  D(z)  is  physically  un¬ 
realizable.  By  increasing  D(z)  becomes  realizable 

although,  obviously,  the  cost  increases.  One  possible  physically 
realizable  solution  can  be  obtained  by  choosing  the  minimal 
possible  S-j  for  which  =  0;  that  is  =  5L^.  Therefore, 

us ing  Theorem  4.7, 
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D(z)  =  -  (54) 

zp  (z)p  Cz)ttCz) 

where  o^Cz)  and  (z)  are  the  solution  pair  of  equation  (37). 

To  find  the  minimum  squared  error,  substitute  equations 
(50)  and  (52)  into  equation  (27)  to  get 


J 


1 

2^j 


X  q'  (z)p'*'(z)Tr(z) 

J  t‘^(z)p  (z)q  (z) 


q  (z'^)p'^(z'^)TT(2  ^)  t'*'  (z)  t'  (z)  ^ 

■’^^(z’^)p' (z'^)q' (z'^)  P’^(z)p'(z)  ^ 


(55) 


where  by  equation  (20) 

q’  (z‘^) 
q'  (z’^) 

Also,  it  is  possible  to  prove  (see  Appendix  A,  Claim  (A-2)) 
that 


6  . 

=  z' 


q’  (z) 


6  - 
=  q 


q*  (z) . 


p"'(z'^)T~  (z)  ^  ^ 

P  ■  ( z )  X  "^  ( z  ■  ^  ) 


and  equation  (54)  thus  becomes 
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Again,  with 


we  have 


1  nCz)nCz  )  dz 

^  p' Cz)P*  Cz'^)  ^ 


“1  TT  ^ 

TT(Z  =  2  TT(Z) 


p'(z'^)  =2  P  (Z) 


j  =  1  ({  .  dz 


J  p‘  (z)p'  (z)z 


Hence , 


J  =  I  RescV"^^^^  - ] 


z  ^p' Cz)p'  (z) 


where  the  summation  is  on  all  the  residues  inside  tlie  unit 

- 

circle;  i.e.  at  the  poles  of  z  p  (z) . 


( 1 )  Consider  the  plant 


G(z) 


and  the  transformed  step  input  sequence 


U(z)  =  . 

From  equation  (1)  we  get, 

q^(z)  =  z  -  0.5,  q  (z)  =  z  -  1,2 

p^(z)  =1,  p  (z)  =  z  -  2 

and  by  equation  (22) 

=  (z-l)a-z)  ' 

Hence , 

x'^Cz)  =  z  x‘(z)  =  1 

p'^Cz)  =  z  -  1  p'  (z)  =  1  -  z. 

Equation  (37)  becomes 

(z-2)a(z)  +  (z-l,2)z(z-l)Tr(z)  =  z  •  (1- 2z)  (1  - 1 , 2z)  • 

where  we  chose  =  1 .  Dividing  both  sides  by 

have,  in  the  notation  of  Lemma  4.4, 
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r 


Cz-2)a'(z)  +  (z-1,2)  Cz-l)7fCz)  =  z  (1- 2z)  (1- 1 , 2z) 

The  unique  solution  of  this  equation  that  satisfies 

=  '^a  '  ^  ^  1  1  and  6^  =  6  _  +  S,  -  1  =  1  is 

T  q  p  ‘ 

given  by 

Ti(z)  =  5.7  +  2.4z 
a' (z)  =  3.42  -  3.62z. 


The  compensator  D(z)  (equation  (54))  is,  therefore, 


nr7l  -  (z-1/2) (5.42-5.62Z) 
(z-1)  (5.7  +  2.4Z) 


and  using  equation  (56),  the  minimum  square  error  is 


J 


I 

z  =  0 


Res  [ 


(5.7-t-2.4z)  (2.4-*-5.  7z)  1 
z(z-2)(i-2z)  ^ 


z 


17.31. 


(2)  Consider  the  plant 


G(z) 


(z-3/2) (z+2) 
(z-a) (z- B) 


l“U  I  3|  <  1 


and  the  transformed  sequence 


U(z) 


Cz-3) 

z-1 


-r  ■  iiiMiMf  imliliifcll  ill  tl 


A 
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For  this  case  we  have, 

q'^Cz;)  =  (z-a)(z-3),  q'Cz)  =  1 

P^(z)  =1  ,  p'  (z)  =  (z-3/2)Cz  +  2) 

T^(z)  =  1  -  3z  ,  T  (z)  =  z  -  3 

p'^Cz)  =  z  -  1  ,  P  (z)  =  1  -  z, 

and  the  equation  to  solve 

(z-3/2)  (z  +  2)ct(z)  +  z(z-l)7r(z)  =  (l-3z)Cl  -  |z)  (l  +  2z) 

with  6^  =  1  and  6^^  =  1 .  The  unique  solution  can  be  found  by 
equating  different  powers  of  z,  yielding 

a(z)  =  -  i  (l+5z) 

TiCz)  =  j  (ll  +  16z). 

With  this,  equation  (54)  yields  the  compensator 


D(z) 


(z-u) (z-B) (l+5z 
2z(z-i) (11+16Z) 


with  minimum  square  error  J  =  27.55. 


VII.  Conclusions 


Output  feedback  with  a  series  compensator  were  considered, 
in  controlling  the  SISO  discrete  system  to  track  a  given  input 
sequence.  MSE  criterion  of  design  was  used  in  a  variational 
problem  whose  solution  yielded  the  required  compensator  trans¬ 
fer  function  D(z) .  This  compensator  was  constrained  to  be 
physically  realizable  and  the  resultant  closed-loop  system  to 
be  both  stable  and  physically  realizable. 

It  turns  out  that  for  D(z)  to  be  the  required  optimal 
solution,  it  must  include  two  polynomials  a(z)  and  6(z) 
that  have  to  satisfy  some  polynomial  equation  whose  parameters 
are  obtained  from  the  input  transform  and  the  plant.  On  the 
other  hand  using  the  solution  of  this  polynomial  equation 
Cwhich,  under  certain  polynomial  degree  conditions  is  unique) 

DCz)  is  the  required  optimal  compensator.  Two  examples,  that 
were  worked  out,  present  the  detailed  steps  required. 
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APPENDIX  A 


The  z 


transform  of  the  input  sequence  is  given  by 


U(z)  =  z 


TT  (z-a.)Tr  Cz-“-k) 

_  L  i=l  ^  k=l  ^ 


N 


M 


Tr(z-b„)Tr(z-3^) 


n=l 


m=l 


(A-1) 


where  and  are  nonzero  scalars,  which  satisfy 


Sfl  5  1.  Ibjil  5  1.  I“kl  ^  ^  ^ 


L  is  some  integer  (positive  or  negative).  For  U(z)  to  be 
causal  we  must  have 


N+M>L+I  +  K. 


(A-2) 


Hence , 


$(z)  A  U(z)U(z'^) 


TT (z-a^  )1T (l-«kZ)Tr (z-“v)Tr (1-a.  z) 

i  =  l  ^  k=l  K=1  ^  i  =  l  ^ 


N 


M 


M 


N 


Tr(z-b„)Tr(i-6„^)-n-(.-6^)Tr(i-b„3) 

n-l  m=i  m=l  n=l 


(A- 3) 


where 


S  =  N  +  M  -  (I  +  K)  . 


(A-4) 
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If  we  separate  the  poles  and  zeros  to  those  inside,  and  those 
outside  the  unit  circle,  i.e.,  if 


(j)  (z)()>  (z) 


then  we  must  have 


^'(z)  = 


P^Z) 


f  [1+sgn  S]  I  K 

=  2  TT (z-a- )Tr (l-'^vz) 

i=l  ^  k=l 

K  I 

=  TT  (z-“u)Tr  (1-a.z) 
k=l  ^  i=l  ^ 

7  [sgn  S-1]  N  M 

-  2^  TT  (z-b.jTT  a-B„2) 

n=l  ^  m=l  ^ 


p'Cz)  =  Tr  (z-6„)Tr  (l-b„z). 


I"  n 

m=l  n=i 


(A-5) 


(A-6) 


lA-7) 


(A-8) 


Notice  that  by  equations  (A-2)  and  (A-4),  S  >  L,  so  that  for 
L  >  0,  S  is  non-negative. 


Claim  A-1: 


P^Cz) 


=  'r'(z) 


=  p’(z) 


(A-9) 


(A-10) 


<5  .  . 


(A-11) 
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Proof.  By  equation  (20) 


^■^(z)  =  z  ^ 


and  from  equation  (A- 5) 


^  =  I  [1+sgn  S]  +  I  +  K. 


Thus,  using  equation  (A- 5),  we  have 


T^(z)  =  z^""^  ]T(z’^-ai)Tr(l-“vz'^)  , 
i  =  l  ^  k=l 


and  comparing  with  equation  (A-6),  equation  (A-9)  follows. 

A  similar  proof  applies  to  equation  (A-10).  Finally, 
substituting  equation  (A-4)  into  equation  (A-12),  we  have 


+  =  I  [1  +  sgn  S]  +  N  +  M  -  S 


=  y  [sgn  S-1]  +  N  +  M 


and  comparing  with  equation  (A- 7),  equation  (A- 12)  follows. 


(A-12) 


Claim  A- 2 ; 


P'*'(z'^)'t'Cz) 
p’  Cz)t'^Cz'^) 


(A-13) 
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Proof.  By  equation  (20) 


-<S 


x'^Cz'^)  =  z 


f^(z) 


and  using  equation  (A-9) 

-<5  + 

x^z'^)  =  z  ^  -x'Cz).  (A-14) 

Similarly,  applying  equation  (20)  to  P^(z  ^)  and  using 
equation  (A-10)  we  have 

+  -1  - 

p  (z  ^)  =  2  -P  (z).  (A-15) 

Using  equations  (A-14)  and  (A-15),  together  with  (A-11),  yield 
our  claim. 
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APPENDIX  B 

Consider  the  space  P  of  pairs 

of  polynomials 

1 

(“(z),3Cz))  of  degrees  6^  and  6^ 

respectively,  i.e. 

( 

. 

\ 

a(z)  =  I  a 
i  =  0 

(B-1) 

1 

t 

S  . 

3(z)  =  I  3.Z 

(B-2) 

‘t 

i 

O 

II 

•H 

(where  and  3^  are  numbers)  such  that 

6=6  +  6.  > 

''  q-  »  - 

P 

(B-3) 

where 

v(z)  =  p'(z)a(z)  +  q*(z)3(z). 


Clearly,  P  is  a  linear  space 

The  first  variation  of  J. 

Let  (hj (z) ,h2 (z) 3  €  P  be  movements  of  (®(z),3(z)). 
The  corresponding  increment  of  the  functional  J  of 
equation  (26)  is, 


tfWUiiitfiMiitfiliAiU 
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Aj  =  j(a+h^,e+h2)  - 


=  ^  ^  (l-T*  (z))(l-T'  (z‘^))nz) 


dz 

z 


2Trj 


(p  Cl-T(z))(l-T(2'^))<I>(z) 


,  dz 


where  T(z)  is  given  by  equation  (10)  and 


T’(z)  = 


p'(z)  (a(z)  +  h^(z)) 


p  (z)(oi(z)  +  h^(z))  +  q'(z)(3(z)  +  h2(z)) 


(B-4) 


Using  Taylor's  theorem  we  get 


J  =  •  2^  f  (l-T(z))p[TCz'l)]<I>(z)^ 

-  ^  (f  Cl-TCz'^))yITCz)JHz)^ 

-  ^  (l-T*(z))y^[T*(z-^)]Hz)^ 

-  ^  (l-T*(z'^))y^T*(z)]$(z)^ 

+  ^  j  y  [T*(z)]y[T*(z'^]$(z)^  +  ...  (B-5) 

In  terms  of  the  increment  polynomials  hj^(z)  and  h2(z),  it 
can  be  shown,  that 


- 
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pT(z)  =  E_(.zi_q'Cz),  y[KCz)], 
[vCz)]*^ 


u2x(z)  =  ijP  U)q 
[v(z)]-' 

,  P~(z..)q;(z) 
[v(z)]" 

where 


v(z)  =  p  (z)a(;z)  +  q  Cz)  B(z)  , 
yv(z)  =  p‘(z)h^(z)  +  q'Cz)h2(z), 
uK(z)  =  0(z)h^(z)  -  ot(z)h2(z), 


and 


y^KCz)  =  0. 

The  (*)  sybolize  the  fact  that  the  corresponding  terms  are 
evaluated  at  i»(z)  +  9j^hj^(z)  and  3Cz)  +  02h2(z)  with 
0  <  0j,  02  <  1.  It  is  easy  to  show  by  substituting  z  ^  for 
z  that  the  first  and  the  second  integrals  of  equation  (B-5) 
are  identical  and  so  are  the  third  and  fourth.  Therefore; 
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^  ^  j  ^  I|lp)  C3(z'^)h^(z-1)  -  a(z-l)h2(z'^))^ 


^  4  f  ^--V^  ^(z)  (p'Cz'Sh  (z'l)  +  q'(z'^)h  (z'b) 

J  J  V  (z)  [v*(2  1  2 


•  (B*(z‘^)hj(z‘b  -  (^*(z'bh2(z'^))^ 

1  I  $(z)  •  P-C^  CB*(z'^)h  (z'b  -  a*(z'^)h  (z'^)) 

^  ^  [v*(z)]^  [v*(z  1  2 


dz 


.(B  (z)h^(z)  -  a  (z)h2(z))^ 


(B-6) 


1 


where 


a*(z)  =  a(z)  +  Bjh^(z) 

6*(z)  =  3(z)  + 

for  some  0  <  9^^,  9^  <  1  and 

v*(z)  =  p‘(z)a*(z)  +  q'(z)B*(z)  • 

By  the  definition  of  the  linear  space  C“*  (z)  ,  B*  (z) )  P. 
Also,  <5  *  =  <5y  and  6[p'(z)h^(z)  +  q'(z)h2(z)]  =  6^.  Using 
equation  (20)  ,  we  have  for  the  first  integral  in  (B-6) 
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[v(z)  ] 


z  ^C6(z)h, (z) 


-  a(z)h2Cz))^ 


(B-7) 


where  we  used  the  following 


26  (6  +6  +  6  +  6  )  =  2(6  +  6  )  -  (6  +6  +  6^  +  6  ) 
Vpqa^  q*"  Pq 


=  6  .  +  6^  -  (6  ,  +  6^). 
q  p 


For  the  second  integral  in  (B-6)  we  have, 


AJ 


2 


q'..CiJ  glCjI  $  (z) 

v*(z) 


p~  (z)q~  (z) 
[v*(z)]^ 


(p‘ (z)K^ (z) 


T  j 

z  (3*(z)h^(z)  -  a*(z)h2t^^)^  • 


q’ (2)^2  ) 

(B-8) 


Using  the  norm 

'a 

IlCh  ,h2)||  =  max(  I  I  h,  •  1  ,  I 
1  ^  i=0  i=0 

AJ2 

it  is  possible  to  show  that,  |[  ^  ||  0  as  l|(hj,h2)|l  0. 

* '  X  2 


For  the  second  integral  in  (B-6)  we  have, 
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1 ,  then 


for  any  P(z)  of  a  given  degree  6^^ ,  <5^  ^  - 

A(z)/B(z)  must  be  analytic  inside  the  contour  of  integration. 


Proof.  By  contradiction,  assume  that  B(z)  has  a  pole  of 
multiplicity  v,  v  <  6g  inside  the  contour  of  integration, 
i.e.,  B(z)  =  (z)  (z-u)'^,  Cj^(ot)  f  0  for  some  ct.  Choose  U(z), 

to  have  a  pole  at  “  of  multiplicity  v  -  1,  i.e., 
u(z)  =  C2(z)(z-a)  C2(“)  0.  Hence,  by  the  residue  theorem 


U (z)dz 


A(a)C2(t^) 


0 


which  contradicts  our  assumption. 
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